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Abstract 

We provide a unified investigation of normal modes and sound propagation 
at zero temperature in Bose condensed gases confined in highly asymmetric 
harmonic traps and interacting with repulsive forces. By using hydrodynamic 
theory for superfluids we obtain explicit analytic results for the dispersion law 
of the low energy discretized modes for both cigar and disk shaped geometries, 
including the regime of large quantum numbers where discrete modes can 
be identified with phonons. The correspondence with sound propagation in 
cylindrical traps and the one-dimensional nature of cigar type configurations 
are explicitly discussed. 
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The recent experimental realization [|IJ of Bose-Einstein condensation in dilute gases of 
alkali atoms confined in magnetic traps has opened unique perspectives in the study of the 
dynamical and statistical behaviour of mesoscopic quantum systems. These highly degen- 
erate quantum gases have been shown to exhibit both the dynamic features of microscopic 
many-body systems, associated with the occurrence of quantized collective excitations (nor- 
mal modes), and the ones of macroscopic quantum fluids, characterized by the propagation 
of sound in the collisionless regime (zero sound). 

The frequencies of the collective excitations have been the object of various experimental 
and theoretical studies, based on both analytic |^|-f| and numerical |J investigations. 
The magnitude of these frequencies is fixed by the oscillator frequency of the trapping 
potential (typically a few ten Hz), the exact value depending on the nature of the excitation 
(angular momentum etc.) and on the effects of two-body interactions. In the experiment 
of ||] the large number of Bose condensed atoms and the elongated geometry of the trap 
have allowed for very precise "in situ" images of the oscillations of the axial radius and 
consequently for high precision in the measurement of the collective frequencies, whose 
values turn out to be in excellent agreement with the predictions of theory ||. These 
experiments, however, do not probe directly the phononic nature of the excitation because 
of the discretization of frequencies imposed by the harmonic trapping. 

The fact that phonons can propagate in the medium in a continuous way, according to 



Bogoliubov theory, has been recently shown in the experiement of jnj where wave packets 
have been produced and directly observed in a Bose condensed gas of sodium atoms confined 
in a highly asymmetric, cigar type trap. 

The purpose of this letter is to discuss the correspondence between these two dynamic 
features (occurrence of discretized collective modes and phonon-like propagation) and to 
give analytic results for the dispersion law of the collective modes in the case of highly 
deformed traps, including the regime of large quantum numbers where the discrete modes 
can be identified with phonons. 

The proper theory to describe the dynamic behavior of interacting Bose gases is the 



time dependent Gross-Pitaevskii |ll[| equation for the order parameter. This equation can 



be written in a convenient form by expressing the order parameter in terms of its modulus 
and phase, $ = y/pe 1 ^, and looking for equations for the density p and for the velocity field 
v = V0. The equations of motion then take the following form ||: 

^p + V(vp) = (1) 

and 

m|v + V(V„ 4- gp - ^VVJi + imv») = (2) 



where g = 47r ^ ° is the interaction coupling constant, fixed by the s-wave scattering length 
a and 

t 7- 1 2 2,-'- 2 2 /qi 

is the trapping potential, for which we have made the choice of an axially symmetric oscillator 
(rj_ = \Jx' 2 + y 2 is the radial coordinate). Eq.(|l]) is the equation of continuity, while @ 
establishes the irrotational nature of the motion. These equations have the typical structure 
of the dynamic equations of superfluids at zero temperature (see for example [|T2|1 ). 

If the repulsive interaction among atoms is enough strong, than the density profiles 
become smooth and one can safely neglect the kinetic pressure term proportional to h 2 in 
the equation for the velocity field. This yields, for the static solution of (0), the so called 
Thomas- Fermi expression 

flj(r) = i(Ai-VUr)) (4) 

if /i > V ex t(r), and zero elsewhere. Here \x is the chemical potential, fixed by the normaliza- 
tion of p(r). In the case of harmonic trapping one has 

p = hu (l5N—y/ 5 (5) 
l duo 

where luq = (lu z lu 2 _) 1 / 3 is the geometrical average of the oscillator frequencies and ano = 



% /mojQ is the corresponding oscillator length. The Thomas- Fermi approximation (|J) for 



the ground state is accurate to the extent that the conditions \x ^> hu z ,huj± are satisfied. 
The density (f|) has the shape of an ellipsoid with radial (R±) and axial (Z) radii defined 
by muj\R\ = mu 2 Z 2 = 2/i. 

In the following we will neglect the kinetic energy pressure term also in the solution of 
the time dependent eqs. ([I], 0). After linearization these equations then take the simplified, 
hydrodynamic form 

^ = v(c 2 (r)Wp) (6) 
where 5p = p(r,t) — po( r ) and c(r), defined by 

mc 2 (r) = fi- V ext (r), (7) 



can be interpreted as a local sound velocity. The validity of the HD equations (|6|,[7]) is based 
on the assumption that the spatial variations of the density are smooth not only in the ground 
state, but also during the oscillation. In a uniform system (V ex t = 0) this is equivalent to 
imposing that the collective frequencies should be much smaller than the chemical potential. 
In this case the solutions of the HD equation (0) are sound waves propagating with the 



Bogoliubov velocity c = J p/m where /i = gp and p is the equilibrium density. Due to the 
non uniform nature of the trapping potential the solutions of (^,0) exhibit new interesting 
features. The corresponding spectrum is discretized and its explicit form, in the case of 
spherical trapping, was first derived in ||. The fact that these equations have analytic 
solutions reflects the occurrence of underlying, non trivial symmetries as recently discussed 
in ||. For the lowest multipolarities it is possible to obtain analytic results also for arbitrary 
values of u z and uj±. In particular, for the lowest m = value of the z-th component of 
angular momentum and even parity, one finds || 



uj 2 ( m = 0) = 2u 2 ± + -u 2 z =F - yj9u* - 16u 2 ujI + 16u{ (8) 

The dispersion (||) has been also derived in using different approaches. 

In the limit of highly deformed, cigar shaped geometry {oj z u>±), Eq.(|8|) gives the 



result u = Jb/2uj z and 2u± for the low and high energy solutions respectively. As antici- 



pated in the introduction, in the experiment of ||], carried out on a very asymmetric trap 
(uj z /uj± = 17/230), it has been possible to measure at low temperature the frequency of 
the low energy m = mode of even parity with very high precision (u) exp = 1.569(4)co> 2 ) 
in excellent agreement with theory = 1.58). This agreement is not a surprise be- 

cause the conditions of applicability of hydrodynamic theory are very well satisfied in this 
experiment. In fact the parameters of the trap (a/ano ~ 10 -3 ) and the number of atoms 
(N ~ 10 7 ) are such that p ~ 30hu± ~ 400huj z and hence the validity of the Thomas-Fermi 
approximation is well ensured. Furthermore, since the lowest frequencies are of order uj z the 
condition Tius -C p., relevant for the applicability of hydrodynamic theory, is also very well 
satisfied. 

It is useful to derive the dispersion law also for the excitations with higher quantum 
numbers in cigar shaped traps. The excitations we are interested in have frequencies of 
order u z , much smaller than u±. Let us rewrite the linear equations (|||7|) in the form 

uj 2 5p = -—V± ((/i - V aBt )V±Sp) - — V 2 ((// - V ext )V z 5p) . (9) 
m m 

This equation shows that in the limit u z u±_ the low energy solutions cannot have any 

dependence on the radial coordinates. This would in fact yield high frequencies components 

of order in the solution, due to the first term in the r.h.s of @. It is then natural to 

expand the m = solutions of (§) in the form 

5p(r)=6p (z) + \ 2 r 2 ± 6 Pl (z) + ... (10) 

where A = lj z /u)± is the deformation parameter of the trap. After inserting (0) into the 
hydrodynamic equation (|) and integrating over the radial coordinates we obtain, for A — > 0, 
the following differential equation for 5po: 

u; 2 6 Po (z) = -\u 2 z {Z 2 - z 2 )V 2 z 6p (z) + u 2 z zV z 6 Po (z). (11) 

where — Z < z < Z. The discretized solutions of ( [TT]) are polynomials of the form 8p^\z) = 
(z k + az h ~ 2 + ..), satisfying the orthogonality condition dziZ 2 — z 2 )5pq\z)5pq \z) for 
k ^ k'. They obey the dispersion relation 
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<J* = -k{k + 3)<J 2 K . (12) 

already derived in M using a different approach. The number of nodes of these solutions is 
equal to k/2 for even k, and to (k + l)/2 for odd k. 

It is also interesting to look for solutions of (|ll|) localized in the center of the trap (z ~ 0). 
These are sound waves propagating with velocity 

C W = vW 2m ( 13 ) 

where we have used the identity p = ^mu 2 Z 2 for the chemical potential. Notice that in 
the Thomas-Fermi approximation the chemical potential is always related to the central 
density by the relation p = gpo(0) (see eq.(f|)), so that the sound velocity c\d is smaller by a 
factor \/2 with respect to the Bogoliubov velocity calculated in the center of the trap. The 
occurrence of this factor was first pointed out in |13[ and has a simple physical meaning. In 
fact, in deriving the relevant hydrodynamic eq.(|TT|), we have integrated (^) over the radial 
variables, so that the new sound velocity corresponds to an average whose value is smaller 
than the one in the center of the trap. 

To better understand the propagation of sound waves in the case of highly elongated 
traps let us consider a trap with cylindrical geometry and harmonic confinement in the 
radial direction. The hydrodynamic equations in this case are simply obtained by setting 
uj z = in (0) and take the form 

u?6p = -iwiVj. ((Rl - rl)VJp) - l -u\{Rl - rl)Vl8p (14) 

defined in the interval —L < z < L, where 2L is the length of the cylinder, and < 
r_i_ < R±. It is worth noting that in the cylindrical geometry the validity of the Thomas- 
Fermi approximation for the ground state is guaranteed by the condition p ^> hu± or, 
equivalently, Na/L ^> 1 (in this case we always assume L 3> R). If we impose periodic 
boundary conditions at z = ±L the solutions of ( |14| ) can be written in the form Sp = 
(Sp (z) + r\5pi(z) + ..) with 5p and 5pi proportional to e iqz . After integration in the radial 
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variables eq.fll4|) takes, to the lowest order in q 2 Rj_, the simplified form u 2 5po = — ^m^lSpo, 
yielding the dispersion u = cmq with the sound velocity given by (|T3|). The numerical 
solution of (|T^) with larger q has been carried out in jnj). It is not difficult to calcualate 
the first correction to the linear behavior. One finds Spi(z) = q 2 5po(z) /8 and 

u 2 = c 2 1D q 2 (l - ±q 2 Rl) (15) 

Result (pL5f) explicitly reveals that the linear dispersion holds if the wavelength is much larger 
than the radial size of the condensate and that the sound has a negative dispersion [13]. 



Coming to the dynamic behavior in the presence of harmonic trapping one expects to 
observe wave packets propagating with the ID sound velocity c\d if the conditions qZ ^> 1 
and qRj_ <C 1 are simulataneously satisified. Of course the condition Tiq <C mc must be also 
satisfied because it ensures the applicability of hydrodynamic theory. The condition qZ ^> 1 
guarantees that the medium can be treated as locally uniform in the z direction and that 
one can consequently observe wave packets propagating in the central region of the trap. 
The condition qR± <C 1 instead ensures that we are not exciting the motion in the radial 
direction and that the dispersion will look "one dimensional" and given by the first term of 
(Pip. In the experiment of [TTJ Z is a few hundred microns, R± = XZ is a few ten microns 
and mc/h ~ 2 — 4(/im) _1 depending on the value of the peak density. It is then possible 
that the wave packets observed in JlO| are, at least partially, built up with wave vectors 



satisfying the above conditions. A detailed discussion of the propagation of wave packets, 



with the inclusion of non linearity effects, has been recently reported in [14[ 

Let us complete our discussion by calculating the first corrections to the dispersion 
relation ([12]). By solving the HD equations ( |i~4"D to the next order in A 2 we obtain, after 
some straightforward algebra, the result Spi(z) = —^V 2 z 5p (z) (notice the analogy with the 
result 5 pi = q 2 Sp /8 holding in cylindrical geometry) and 

u 2 = h(k + 3K 2 (l - £{k -l)(k + 4)) . (16) 

Some remarks are in order here. First one recovers the limit of the sound wave dispersion 



(jnj|) holding in cylindrical geometry in the limit of large quantum numbers k 3> 1, by the 
proper identification 

cW = kW z /A (17) 

yielding k 2 = q 2 Z 2 . This is consistent with the already discussed condition qZ S> 1 needed 
to observe phonons propagating in the z-th direction. In the same limit the first corrections 
in the dispersion and (|16|) coincide since one has q 2 R 2 j_ = k 2 X 2 . This completes the 
correspondence between the propagation of discretized modes and sound and shows the 
analogy between the dynamic behavior in the cylindrical and elongated harmonic oscillator 
geometries. 

Concerning the frequencies of the discretized modes predicetd by (|1^) it is worth pointing 
out that the lowest mode (k = 1) corresponds to the center of mass motion and its frequency 
coincides with the oscillator frequency u z . This frequency is unaffected by the presence of 
two-body interactions. The second mode (k = 2) is the "quadrupole" collective excitation 
observed in [|J. It corresponds to the low energy solution ([§]) in the A — > limit. The 
direct experimental observation of the higher modes, as well as of the first correction in A 2 
predicted by ( |T6"D would complete the scenario of the low energy excitations in the elongated 
geometry. 

It is finally interesting to discuss the one-dimensional nature of these systems. One 
should first point out that all the results discussed in this paper have been derived starting 
from 3D configurations. In particular, in order to derive the dispersion law ([16]), we have 
assumed the validity of the Thomas Fermi approximation in both axial and radial directions, 
so that (|16"D holds if /i 3> hu± 3> Tiu z . This means that the ground state wave function of 
the system is built up including many excited single particle states in both axial and radial 
directions. A full ID problem should involve only the lowest oscillator wave function in the 
radial direction and in this case the corresponding excitation spectrum in the Thomas- Fermi 
regime would be |T5[ u 2 = jk(k + l)u 2 instead of ([11]). This dispersion is easily derived from 



(|) ignoring the radial coordinates in the equation and holds if fkv± ^> (/x — hu±) ^> fao z . 



At present this regime is far from experimental possibilities. Nevertheless, even remaining 
in the 3D Thomas-Fermi regime, it is clear that for highly elongated configurations the low 
energy dynamic behavior (qR± <C 1) looks one-dimensional, the radial directions provid- 
ing only a renormalization of the sound velocity. So all the statistical and thermodynamic 
properties of ID systems should apply to these configurations provided the temperature is 
smaller than the radial oscillator energy. This includes in particular the Luttinger liquid-like 
behavior recently suggested for these systems |16| and the two-step Bose-Einstein conden- 



sation recently discussed in the context of the ideal Bose gas [fL7|1 . Furthermore, due to the 
coupling with the radial modes, the transition from the phonon to the single particle regime 
exhibits new interesting features. In particular the first correction to the phonon disper- 
sion is negative (see (|15|) and (0)), differently from the traditional Bogoliubov behavior 
uj 2 = c 2 q 2 + (q 2 /2m) 2 . 

In analogous way we can carry out the analysis in the disk geometry [lo z ^> uo±). In 
this case, to the lowest order in 1/A 2 , the density fluctuations will depend only on the radial 
coordinates and, after integration of Eq.(|S]) in the variable z, the relevant equation for 5p(r^) 
takes the form 

uJ 2 5p(r ± ) = -^l((R 2 ± - r 2 ± )V 2 ± 6p(r ± ) + u 2 ± v ± ■ V ± 6p(r ± ) (18) 

Notice that in this case wave packets in the center of the trap will propagate with the " 2D" 
sound velocity 



02D = V|-- (19) 

V 3 m 

The discretized solutions of (|T8"D) have the form 

Sp (n,m) = ( r 2n + + __yrnjm4 ( 2Q ) 

where m is the z-th component of angular momentum and n fixes the number of radial 
nodes. They satisfy the orthogonality condition J r±< R ± dr±(R 2 L — r^) 1 / 2 ^ p("' m ) '5 p^ n ' ,m '^ = 
for n,m ^ n', m'. The resulting dispersion takes the form 



uj 2 = {-n 2 + -nm + 2n + m)uj_. (21) 



The cases n = 0,m = 2 (a; = 2c<j_i_), and n = l,m = (u = yl0/3cjj_) corespond to the 
modes of even parity observed in the experiments of (the m = state corresponds to 
the low energy solution (|8|) in the u)± uj z limit). One should however note that in this 
experiment the deformation of the trap and the number of atoms in the condensate were 
not very large (A = y/8 and iV ~ 10 4 ). As a consequence the conditions required to apply 
the dispersion law ( f2~l~D (A > 1 and validity of the Thomas-Fermi approximation) are not 
very well satisfied in this case. 



We finally note that also in deriving the dispersion law (21) we have assumed the validity 
of the Thomas- Fermi approximation in both axial and radial directions (/x 2> huj z 3> hu±). 
The hydrodynamic dispersion in a true 2D trap would in fact follow a different dispersion 
law given by [|TjJ uo 2 = (2n 2 + 2nm + 2n + m)uj_. This 2D hydrodynamic dispersion law 
holds if the conditions hu z 3> (/i — \Tuo z ) tuv± are satisfied. 

Stimulating discussions with Dan M. Stamper-Kurn and Wolfgang Ketterle are acknow- 
eledged. 
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